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We develop a theory of cooperative emission of light by an ensemble of emitters, such as fluo- 
rescing molecules or semiconductor quantum dots, located near a metal nanostructure supporting 
surface plasmon. The primary mechanism of cooperative emission in such systems is resonant energy 
transfer between emitters and plasmons rather than the Dicke radiative coupling between emitters. 
We identify two types of plasmonic coupling between the emitters, (i) plasmon-enhanced radiative 
coupling and (ii) plasmon-assisted nonradiative energy transfer, the competition between them gov- 
erning the structure of system eigenstates. Specifically, when emitters are removed by more than 
several nm from the metal surface, the emission is dominated by three superradiant states with the 
same quantum yield as a single emitter, resulting in a drastic reduction of ensemble radiated en- 
ergy, while at smaller distances cooperative behavior is destroyed by nonradiative transitions. The 
crossover between two regimes can be observed in distance dependence of ensemble quantum effi- 
ciency. Our numerical calculations incorporating direct and plasmon-assisted interactions between 
the emitters indicate that they do not destroy the plasmonic Dicke effect. 



I. INTRODUCTION 

Superradiance of an ensemble of dipoles confined 
within a limited region in space has been discovered in 
the pioneering work by Dicke^i The underlying physical 
mechanism can be described as follows. Suppose that a 
large number, TV, of dipoles with frequency wq are con- 
fined in a volume with characteristic size L much smaller 
than the radiation wavelength Aq = 2t: /loq. Then radia- 
tion of an ensemble is a cooperative process in which the 
emission of a photon is accompanied by virtual photon 
exchange between individual emitters. This near field ra- 
diative coupling between the dipoles leads to formation of 
new system eigenstates, each comprised of all individual 
dipoles. The eigenstates with angular momentum / = 1 
are superradiant, i.e., their radiative lifetimes are very 
short, ~ t/N, where r is radiative lifetime of an indi- 
vidual dipole; the remaining states are suhradiant with 
much longer decay times, ^ t{Xo/L)'^ 3> t. 

Since the appearance of Dicke paper, cooperative ef- 
fects based on Dicke radiative coupling mechanism have 
been extensively studied in atomic and semiconductor 
systems (see, e.g., reviews in Refs. [2-13) • Two differ- 
ent decay times corresponding to superradiant and sub- 
radiant states were observed in a system of two laser- 
trapped ionsS- and, more recently, in laterally arranged 
quantum dotsi^ Other examples of cooperative behavior 
analogous to the Dicke effect include, e.g., electron tun- 
neling through a system of quantum dots^*^ and sponta- 
neous phonon emission by coupled quantum dots4»2. 

Recently, we extended the Dicke effect to plasmonic 
systems comprised of N dipoles located in the vicinity 
of a metal nanostructure, e.g., metal nanoparticle (NP), 
supporting localized surface plasmon (SP).^*^ In such sys- 
tems, the dominant coupling mechanism between dipoles 
is plasmonic rather than radiative, i.e., it is based on 
virtual plasmon exchange (see Fig. [Ij. This plasmonic 
coupling leads to formation of collective states, simi- 
lar to Dicke superradiant states, which dominate pho- 
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FIG. 1: (Color online) Radiative coupling of emitters in free 
space (a), and plasmonic coupling of emitters near a metal 
nanoparticle (b). 



ton emission. Furthermore, the nanostructure acts as a 
hub that couples nearby and remote dipoles with about 
equal strength and hence provides a more efficient hy- 
bridization of dipoles compared to radiative coupling. 
In general, as dipoles orientations in space are non- 
uniform, there are three superradiant states with total 
angular momentum 1 = 1, each having radiative decay 
rate ~ where F'' is radiative decay rate of a 

single dipole near a nanostructure (i.e., with plasmon 
enhancement ) 

The principal difference between plasmonic and usual 
(photonic) Dicke effects stems from non-radiative energy 
transfer between the dipoles and the nanostructure. Let 
us first outline its role for the case of a single dipole 
near metal NP. When an excited emitter is located close 
to metal surface, its energy can be transferred to op- 
tically inactive excitations in the metal and eventually 
dissipated (Ohmic losses). This is described by the non- 
radiative decay rate, F"'' cx d~^, where d is the dipole- 
surface separation.— Note that very close to metal sur- 
face (^ 1 nm), this dependence changes to cx d^^ due 
to surface-assisted generation of electron-hole pairs out 
of the Fermi sea.— As a result, the radiation of a cou- 
pled dipole-NP system is governed by a competition be- 
tween non-radiative losses and plasmon enhancementii^ 
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that determines system quantum efBciency, Q = r'"/r, 
where F = r*" + r"*" is the full decay rate. Indeed, the 
radiated energy is = {hkc/2)Q, k and c being wave 
vector and speed of hght, and its distance dependence fol- 
lows that of Q. Namely, with decreasing d, the emission 
first increases due to plasmon enhancement, and then, at 
several nm from metal surface, it is quenched due to sup- 
pression of Q by non-radiative losses. Both enhancement 
and quenching were observed in recent experiments on 
fluorescing molecules attached to a metal NP^i^"— and, 
not too close to NP surface, the distance dependence of 
single-molecule fluorescenceiiii^ was found in excellent 
agreement with single dipole-NP models»^"~ 

When radiation takes place from an ensemble of emit- 
ters near a metal nanostructure, there are two distinct 
types of plasmon-induced couplings between the emit- 
ters. The first is plasmon- enhanced radiative coupling, 
described by radiative decay matrix F^^, where indexes 
j,k — 1, . . . ,N refer to emitters, that is a straightforward 
extension of Dicke radiative coupling obtained by incor- 
porating SP local field into the common radiation field. 
Correspondingly, the eigenstates of F^^ are superradiant 
and subradiant states characterized by the strength of 
their coupling to radiation field. In the ideal case of 
"point sample," i.e., kL <^ 1, the subradiant decay rates 
are negligibly small and FJ^ essentially has just three 
non-zero eigenvalues, corresponding to superradiant de- 
cay rates, each scaling with as ^ NT^ 

The second coupling mechanism is non-radiative en- 
ergy transfer between dipoles that takes place in two 
steps: an excited dipole first transfers its energy to plas- 
mons in nanostructure via its electric field, and then this 
energy is transferred to another dipole. This process in- 
volves plasmons with all angular momenta I, and it is 
described by non-radiative decay matrix, F"^. Impor- 
tantly, plasmons with I > 1 couple to both superradiant 
and subradiant states, so that the eigenstates of full decay 
matrix, Tjk = F^^, -l-F"^, are not superradiant and subra- 
diant states, but their admixtures. Close to metal surface 
where non-radiative processes are dominant, F"^ prevails 
over Tjf, , and no cooperative behavior is expected. How- 
ever, when dipoles are removed from the surface by more 
than several nm, the energy transfer occurs primarily via 
optically active dipole surface SP and therefore no signif- 
icant mixing of superradiant and subradiant states takes 
place and superradiance is intact. 

This observation was confirmed by numerical calcula- 
tion of eigenvalues of Tjk, i.e., full decay rates of sys- 
tem eigenstates, for ensemble of N dipoles randomly 
distributed in a solid angle around a NP.iS Namely, in 
a wide range of dipole-NP distances, three eigenvalues 
corresponding to superradiant states are well separated 
from the rest and scale with N according to ^ A^F/3. 
Since the superradiant states are the optically active ones 
with radiative decay rate also scaling as ~ NT^ /3, their 
quantum efficiencies essentially coincide with Q of single 
dipole-NP system. Therefore, in the cooperative regime, 
the ensemble quantum efficiency, Qens, is thrice that of 



the single dipole-NP system, 

Qens — (1) 

regardless of the ensemble size. Thus, the total radiated 
energy of the ensemble. Wens = {fikc/2)Qens, is reduced 
to just 3W. The remaining energy is trapped by iV — 3 
subradiant states and eventually dissipated in the metal 
rather than being emitted with a much slower rate, as 
it would be the case in free space. On the other hand, 
at several nm from metal surface, the non-radiative cou- 
pling is dominated by higher I plasmons causing strong 
mixing of superradiant and subradiant states, i.e., all 
system eigenstates have comparable quantum efficien- 
cies and Qens oc N. Therefore, with increasing distance, 
Qens, should first exhibit a sharp rise with its slope oc N, 
and then switch to a more slower SQ dependence. 

Indications of such behavior were reported in the re- 
cent experiment by Dulkeith et al.^ where a systematic 
study of the ensemble fluorescence vs. distance to metal 
surface was performed for Cy5 fluorophores attached to 
Au NP in water. The distance was controlled by vary- 
ing fluorophores concentration; with increasing concen- 
tration, the linker molecules stretched outwards to ac- 
commodate repulsive dipole-dipole interaction between 
fluorophores. Therefore, within isolated dipole-NP pic- 
ture, one would expect that for larger distances, at which 
fluorophores concentration was higher than average, the 
measured quantum efficiency, normalized to some aver- 
age fluorophores number, should have exceeded the cal- 
culated efficiency in single dipole-NP modelsi.22, Instead, 
with increasing distance, the normalized Q, exhibited 
rapid saturation and, at large distances, was consider- 
ably smaller in magnitude than the calculated one. 

In the above discussion, we completely ignored inter- 
actions between the dipoles. In fact, the role of dipole- 
dipole interactions in cooperative emission is highly non- 
trivial since they introduce a disorder into the energy 
spectrum by causing frequency shifts among randomly 
distributed in space but otherwise identical emitters.^— 
Mesoscopic cooperative emission from a disordered sys- 
tem, i.e., for particular disorder realizations rather than 
its averaged effect, was considered in Ref. [2^ (see also 
Ref. 0). It was found that frequency shifts due to 
dipole-dipole interactions lift the degeneracy of subradi- 
ant states without having signiflcant impact on superra- 
diant states. It was also shown that interactions between 
collective eigenstates are much weaker than those be- 
tween individual dipoles, including typical nearest neigh- 
bors (i.e., separated by ~ LN~^/^), due to cancellations 
between dipole-dipole terms among individual pairs with 
their constituents belonging to different eigenstatesi^ On 
the other hand, for completely random distribution and 
high concentration of dipoles, the rare instances of ex- 
tremely close dipoles (i.e., with separation ^ LiV~^/^) 
can prevent the formation of superradiant states. 

The role of interactions in plasmon-niediated coop- 
erative emission is characterized by several distinctive 
features. First, since the decay matrices Tjk contain 
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plasmon pole, the relative strength of dipole-dipole in- 
teractions is effectively reduced as compared to purely 
photonic case. Second, there are additional corrections 
to emitters' frequencies, one originating from plasmon- 
enhanced radiative coupling and another from nonradia- 
tive coupling. While neither of those corrections diverges 
as dipoles approach each other, the latter becomes very 
large as dipoles approach the metal surface due to contri- 
bution of high-Z plasmons. Therefore, the actual system 
eigenstates are determined by both the interactions and 
the energy exchange, and their frequencies and decays 
rates must be found simultaneously. It is precisely the 
goal of this paper to calculate the full spectrum of inter- 
acting emitters near a metal NP. 

Specifically, we consider a common situation when 
emitters, e.g., fluorescent molecules or quantum dots, are 
attached to NP surface via flexible linkers. Typically, flu- 
orophores bound to linker molecules have certain orienta- 
tion of their dipolc moments with respect to NP surface 
and, due to repulsive interactions, their angular positions 
are ordered rather then random.— Therefore, we assume 
here that angular positions of emitters coincide with the 
sites of spherical lattice, such as fullerenes. Specifically, 
we perform our numerical simulations for C20, C60, and 
C80 configurations for respective number N of dipoles; 
we also study the effect of deviations from ideal lattice. 
We find that not too close to metal surface, the system 
eigenstates fall into three groups, each dominated by a 
particular coupling mechanism: three superradiant states 
dominated by plasmon-enhanced radiative coupling, one 
state dominated by direct dipole-dipole interactions, and 
the rest dominated by non-radiative coupling via NP. 
Importantly, superradiant states are not significantly af- 
fected by dipole-dipole interactions whose main effect is a 
large frequency shift of subradiant state with the smallest 
decay rate. 

We also address the effect of individual emitters' in- 
ternal nonradiative processes on the ensemble quantum 
efficiency. Internal relaxation is known to inhibit the 
photonic Dicke effect if the corresponding decay rate, 
Fg'', is sufficiently high, i.e., relaxation time is shorter 
than the lifetime of superradiant state<^ However, near 
metal nanostructure, the plasmon-enhanced radiative de- 
cay rate, F'', is significantly larger than Fq'', so that 
plasmon-mediated superradiance is less sensitive to quan- 
tum yield of individual emitters. Specifically, we perform 
numerical calculations for high-yield and low- yield emit- 
ters to show that Eq. ([T|) holds for both types. More 
precisely, in cooperative regime, it holds almost exactly 
for high- yield emitters while for low-yield emitters Qens 
is somewhat larger than 3Q due to the effective iV-fold 
suppression of Fg'' by superradiant states' decay rate. 

An obvious application of plasmon-mediated coopera- 
tive emission is related to fluorescence of a large but un- 
certain number of molecules at some average distances 
from metal nanostructure. For single-molecule case, flu- 
orescence intensity variation with distance was proposed 
to serve as nanoscopic ruler, owing to the excellent 



agreement of measured distance dependences with single- 
dipole models. ^^"^"^ In the case of molecular layer, the am- 
biguities caused by uncertain molecules number and their 
separation from the metal surface prevent, in general, de- 
termination of system characteristics from fluorescence 
variations. However, in cooperative regime, the ambigu- 
ity related to molecules number is removed, and fluores- 
cence intensity is essentially determined by Eq. ([T|) with 
distance-averaged single-molecule quantum efficiency. 

The paper is organized as follows. In Sec. [Hi the 
derivation of plasmonic coupling for an ensemble of 
dipoles distributed near metal NP is given. In Sec. IIII[ 
the general expression for radiated energy is obtained. 
Our numerical results and discussion are presented in 
Sec. IIVI Conclusions and some technical details are pro- 
vided, respectively, in Sec. IVland in the Appendix. 

II. PLASMONIC COUPLING OF RADIATING 
DIPOLES 

We consider a system of N emiters, such as fluorescing 
molecules or quantum dots, with dipole moments dj = 
djej, where dj and are magnitude and orientation, re- 
spectively, located at positions near a metal NP with 
radius R. Throughout the paper, we assume that charac- 
teristic size of the system (NP-|-dipoles) is much smaller 
than the radiation wavelength, |rj— rfe| = \rjk\ <C Aq. We 
also assume that emission events by individual molecules 
are uncorrelated, i.e., after excitation each molecule re- 
laxes through its own internal nonradiative transitions 
before emitting a photon. Then the ensemble emission 
can be described within classical approach by consid- 
ering dipoles as identical Lorentz oscillators with ran- 
dom initial phases driven by common electric fleld, i.e., 
one created by all dipoles in the presence of metal NP. 
The frequency-dependent electric field, E(r,a;), satisfies 
Maxwell's equation 

E(r,c^)-VxVxE(r,c^) = ^j(r,c^), (2) 

where dielectric permittivity e(r, uj) is that of the metal 
inside NP, e(aj) for r < i? and that of outside dielectric eg 

oo 

for r > R. Here j(r,a;) = —i J e^'^*j{t)dt, is the Laplace 



transform of dipole current 

j(t)-g^d,(t)e,<5(r-r,), (3) 
j 

where dipole displacements dj (t) are driven by the com- 
mon electric field at dipoles' positions 

d;+c^2^, = ^E(r,,i)-e, (4) 

with the initial conditions (at < = 0): = doBj simpj, 
dj ~ ujQdoej cosifj, and E = (dot stands for time- 
derivative). Hereafter, wq, g, m, and ipj are oscillators' 
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excitation frequency, charge, mass, and initial phase, re- 
spectively (wq = h/mdQ). Closed equations for dj{(jj) 
can be obtained by using Laplace transform of Eq. dU 
with the above initial conditions and eliminating E from 
Eqs. ([U and (|4]). Laplace transform of Eq. ^ has the 
form 



e{Y,u))(jj 



E(r,w) - V X V X E(r,a;) 

j) [iwodoGj cos (ySj 
j 

—uy^dj {io)ej + ajdoBj sin tpj~\ 



(5) 



At this point, it is convenient to introduce normalized 
displacements 

Vj{u>) = dj{uj)/do — i (ujo/uj'^'j cosi^j — sin^pj (6) 

and the solution of Eq. ^ reads 

'indoquj'^ 

■ 



E(r,L^) 



■^G(r,rj,tj)-ejWj, 



(7) 



where G(r,r',aj) is the electric field Green diadic in the 
presence of NP. From Eq. ([4]), for photon frequency close 
to dipoles frequency, w « wqj we obtain a coupled system 
of equations for normalized displacements, 



E 

fc 



Vk = ^e 



(8) 



where Ejfc = A^fc 
trix, given by 



2^ jk is the complex self-energy ma- 



(9) 



The system Eq. ([5]) determines eigenstates of the ensem- 
ble of N dipoles in common radiation and plasmon field, 
while complex eigenvalues of the self-energy matrix give 
eigenstates frequency shifts with respect to ujq and their 
decay rates. In the absence of NP, real and imaginary 
parts of the self-energy matrix are dipole-dipole inter- 
action and radiation coupling between dipoles j and fc, 
given by (in lowest order in krjk) 



4(fcrjfc) 

Tj J. = Pq ■ Gfc , 

where 
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[ej ■ ek) 



3(ej ■ rjfc)(efc ■ r^-fc) 



3mc^ 3fteo ' 



(10) 



(11) 



is the radiative decay rate of a dipole in a dielectric 
medium, fi = qdo is the dipole moment, and k — y^eow/c 
is the wave vector. The eigenstates of photonic decay 
matrix Pj^j, are superradiant and subradiant states, the 



former having decay rate '-^ -^^ro- In the case when all 
dipoles are aligned, there is only one superradiant state 
that couples to radiation field, while for general dipole 
orientations there are three such states with angular mo- 
mentum ^ = 1. Note that in the longwave approximation 
used here, the decay rates of subradiant states vanish 
("point sample"). 

In the presence of metal nanostructure, the system 
eigenstates are determined by the full Green diadic in 
self-energy matrix Eq. In the case of spherical NP, 
the longwave approximation for G{rj,r^.;uj) can be eas- 
ily foundi^i^'2^ The details are given in the Appendix, 
and the result reads 



jk 



(12) 



,1^ 
2^0 



Jk 



where 



a/(w) 



eiu) + (1 + l/l)eo 



(13) 



is NP ^-pole polarizability. The matrices K^) and t'^^ 



are defined as 



K 



T 



(0 



Ait 



21 



7n— — l 



An 



21 



•-jk 

where 



I SI 

rn— — l 



ij •t/;,™(r,)][efe •x?™(rfc)], (14) 
J • V'im(rj)][efe •i/>;„(rfc)], (15) 

, ximir) = V[r'YUr)], (16) 



and Yimii") are spherical harmonics. For I = 1, these 
matrices can be evaluated as 

J - ek) - 3(ej • fj)(fj - ek) , (17) 
{ej ■ Ok) - 3{ej ■ rj){ek ■ Tj) (18) 

-3(ej • ffc)(efe • ffe) -I- 9(efe • rk){ej ■ rj)(fj • rfe) 
where we have used identities 

1 3 
E '/'i™('^fc)^im(^") = 7-^ - 3ffe(r • rfc)] , 





= — [(e 


J-Jk 


1 r 


-,3,3 [ 



m— — 1 
1 



(19) 



For dipoles oriented normally with respect to NP surface 
(gj = Yj), we obtain 
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(1) 
jk 



,3,3 



(e, • efe), (20) 
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and for parallel orientation (gj • f j = 0) , we similarly get 
-3 (ej • Gfc), T^j,^ — 



(1) 



jk 



(21) 



The decay matrix, Tjk = — 2ImI]jfe, can be decomposed 
into radiative and nonradiative terms, Tjk = F^j, + F"^, 



as follows 



F^^,=r5[(e,.e,)-a;[i^« 

pnr _ 3Fo ^ ^//^(O 



2^(1) 



(22) 



The diagonal elements T^j and F"J describe, respectively, 
plasmon-enhanced radiative decay rate of an isolated 
dipole near a NP, and nonradiative transfer of its en- 
ergy to electronic excitations in metal. The non-diagonal 
elements of F^^ describe plasmon-enhanced radiative cou- 
pling that generalizes the Dicke mechanism responsible 
for cooperative emission by incorporating local field en- 
hancement into near field radiative coupling. On the 
other hand, the non-diagonal terms in F"^ describe non- 
radiative energy transfer between dipoles mediated by 
NP plasmons. The latter coupling is absent in the pho- 
tonic Dicke effect but it plays important role in the plas- 
monic Dicke effect, as we will see below. 

Since the numerical calculations below are carried for 
normal dipoles orientations, we provide here the corre- 
sponding expressions for self-energy matrix T,jk- The de- 
cay matrix has the form 



2"'i (4 + 3- 



4|aip 
^3^3 



2fc3 ^ 



3 ' fe / 
1)2 

2-^P,(cOS7jfe), 



(23) 



where Pi is Legendre polynomial and jjk is the angle 
between positions of dipoles j and k measured from NP 
center (here cos7jfc = ej-ek). The real part of self-energy 
matrix has the form 



A 



jk 







"l — - 









— y 

4fc3 ^ 



^k 



■P;(cOS7jfc) 



(24) 



where the second term describes NP-induced interac- 
tions. The latter in turn consists of two terms, first 
coming from plasmon-enhanced radiative coupling and 
second coming from non-radiative coupling. 

Note that both NP-induced terms are weaker than 
their counterparts in Tjk while having same symmetry 
and therefore is not expected to significantly alter the 
eigenstates. On the other hand, the dipole-dipole inter- 
action term, Eq. (jlOl) . has different symmetry and can 
become large for two dipoles in a close proximity to each 
other. The effect of interactions on cooperative emission 
is studied in Sec. IV. 



III. RADIATED ENERGY 

In this section, we derive general expression for total 
energy radiated by an ensemble of dipoles near a metal 
nanostructure. The radiated energy in the unit frequency 
interval is obtained by integrating the far field (r — >■ 00) 
spectral intensity over solid angle^S 



dW{u!) 
duj 



47r2 



j |E(r,c.)|^2dJ7, 



(25) 



and then averaging the result over initial random phases 
of individual dipoles, Lpj. The electric field E(r.cj) is 
given by Eq. (O, where Vj is the solution of Eq. (jH). 
Then the energy density takes the form 



dW _4r2eoAi^Wo 
duj c-^ 



(26) 



X dnvjvl[G{Y,Yj) ■ e.j\ ■ [G*(r,rfe) • e^] , 

jk 

with G(r, Tj, oj) = G"(r, rj, w) + ^(r, r^, w) replaced by 
its far field asyniptotics (see Appendix), 



Gi.Ar.rj)^ — 



3 



-^f,yi,„(f)xr:„(r, 



m 

where first two terms come from the free space part, G*^, 
and the last term comes from the scattered part, G**, of 
the Green diadic. The angular integral in Eq. ([^5]) can 
be performed using the relations Eq. ([TS]) . The free space 
contribution yields 

4^r2y"dl7[G°(r,r,)-e,] • [G"*(r, r^) • e^] =^(e, -e^). 

(28) 

The other integrals in the product Eq. (^5)) are evaluated 
using the relations 



and 



( — j 2^ ['/'im(i"i) • ej] [i/^i™' (rfc) • Gfe] 

mm' 

X / ^vn™(f).vr;„,(f) 

m 

yielding 

j ^[G(r,r,)-e,].[G*(r,r,).e,] 



(29) 



(1) 



-T 

^ jk ' 



= -K 



(1) 



3 jk ' 
(30) 



(31) 



1 



(47rr)2 3 L 



jk 
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The energy density then takes the form 

2^ Vj^k^jk, 



du! 



(32) 



jk 



where Ajk ~ T, 



\eyek)-a,Kf^-alK^y 



\^i?T^'k' 



Matrix Ajk is not symmetrical, however only its symmet- 
rical part, equal to F^^, contributes to the final expres- 
sion. The solution of Eq. ([8|) can be presented as 



k 



1 



Wo 



(33) 



jk 



and after averaging out over the initial random phases 
Lpj^ we finally obtain 



dW{uj) hkc 



duj 



47r 



Tr 



1 



f 



a; — cjQ ^ St uj — ujo — Y, 



(34) 



where trace is taken over the indexes (jk). The analysis 
of this expression and the results of numerical calcula- 
tions will be presented in the next section. 



S|J) = (Aj — |rj)| J), where Aj is frequency shift of 
collective eigenstate | J) relative to and Tj is its decay 
rate. The molecular relaxation can be accounted for by 
adding to Sjfe a diagonal term, —^SjkTg'^. Then, after 
frequency integration of Eq. (l34l) , the ensemble radiated 
energy takes the form 



ens — „ ^ens, Wens / ^ p T^nr ' 



(37) 



where Fj — {J\t^\J) is radiative decay rate of state | J). 

In the photonic Dicke effect, superradiant and subra- 
diant states are eigenstates of the radiative decay matrix 
obtained from the free space Green diadic. Similarly, 
in the plasmonic Dicke effect, superradiant states are 
eigenstates of plasmon- enhanced radiative decay matrix 
Tjf. . Let us illustrate the emergence of plasmon-mediated 
superradiance for a simple case when all dipoles are at the 
same distance from NP surface and are oriented normal 
or parallel to it. Then it is easy to see that the corre- 
sponding decay matrix, F^^, = ^^^j ■ Gfc with F'' given by 
Eq. ([55]) . has just three nonzero eigenvalues. Indeed, let 
us introduce new decay matrices as 



IV. DISCUSSION AND NUMERICAL RESULTS 



iVF'' 



-B 



(38) 



Let us start with a single dipole located at tq near a 
metal NP. In this case, the self-energy is a complex num- 
ber, S = A - |F, where A = A^-^- and F = r''-hF"'^ = F^.,- 
are single-dipole energy shift and decay rate, respectively. 
For normal (s =-L) and parallel (s =||) dipole orienta- 
tions, using Eqs. (1^ and (1^ . these are given by^S 



where 



F'^ ==F5 



3 



E 



= H1 V 

2fc3 ^ 

^2/-l-4 I ' 



5iV 

2/+4 



(35) 



witha_L = 2, ^ {l + lf andaii = -1, = l{l + l)/2. 
Note that both terms in A are smaller than their coun- 
terparts F'' and r"*" due to plasmon pole in the imaginary 
part of NP polarizability a". Radiated energy of single 
dipole-NP system, obtained by frequency integration of 
Eq. ([M)) . is given by 



W 



hkc hkc 



(36) 



where we included internal molecular relaxation rate, 
Fg'", into quantum efficiency Q. For N uncoupled dipoles, 
i.e., for purely diagonal T,jk — Sjk (A — |F), Eq. ([M)) 
decouples into sum of N independent terms, yielding 

Wens = NW. 

In the presence of inter-dipole coupling, the system 
eigenstates, |J), are those of the self-energy matrix, 
Eq. (jl2p . The corresponding eigenvalues are complex. 



(39) 



is 3 X 3 matrix in coordinate space with Tr_B = 3. It is 
easy to see that Tr[(r'")"'] = Tr[(7'')"] for any integer 
n, i.e., the N x N matrix F^^. has only three non-zero 
eigenvalues coinciding with those of matrix 7^^^ 



iVF'^ 



(40) 



where ^ 1 are eigenvalues of B^^. Note that the decay 
rates of the remaining TV — 3 subradiant states vanish 
in the long wave approximation used here; they acquire 
finite values in the next order in (fcro)^. 

Let us turn to non-radiative coupling, described by 
matrix F"^. Its diagonal elements, F"J, describe nonra- 
diative energy exchange between excited dipole and NP 
plasmon modes with all angular momenta, as indicated 
by polarizabilities a" in Eq. ((23)) . The non-diagonal el- 
ements of F"^ describe a process by which a plasmon 
nonradiatively excited in the NP by dipole k transfers its 
energy to another dipole j. In general, due to high-? plas- 
mons involved in nonradiative coupling, the eigenstates 
of F"^' are different from those of plasmon-enhanced ra- 
diative coupling FJ^ which contains only dipole {I = 1) 
plasmon mode. Therefore, the eigenstates of full decay 
matrix Tjk ~ F"^ + F"^ are not pure superradiant and 
subradiant states but their admixtures. However, the 
high-? plasmon contribution to F"^ is significant only at 
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very small d [see, e.g., Eq. ([23)) ] while for d larger than 
several nm F"^ is dominated by the I = 1 term. In fact, in 
a wide range d, nonradiative coupling between dipoles is 
mainly through the optically active dipole plasmon mode 
that does not cause mixing between superradiant and 
subradiant states. Namely, it can be easily seen from 
Eq. 



that the eigenstates of the / = 1 term in F"^ are 
the same as those of F^j, so the corresponding eigenval- 
ues are similarly given by F^'' — (A^F"'"/3) A^. Thus, the 
superradiant quantum efficiency 



+ FJ,- F + 3FJf77VA^ 



(41) 



only weakly depends on N. Therefore, the sum in 
Eq. ((57)) includes just three terms, yielding 



We. 



hkc 



hkc \ ^ 
~2~ 



^ F- 



SF^VTVA^ 



(42) 



For high-yield (small Fq*") emitters, we obtain Eq. ([T]) 
and hence Wens ~ 3W. In contrast, the radiated power, 



p — 



hkc 



N 



( hkc 



QV = NP, (43) 



scales with the ensemble size due to shorter (by factor 
iV/3) radiative lifetime of superradiant states. 

For low-yield emitters (large Fq""), the relation Eq. ([T]) 
holds only approximately. However, it is evident from 
comparison of Eqs. p6p and ()42l) that here the relative 
effect of internal relaxation is much weaker than for usual 
cooperative emission. Numerical results for both high- 
yield and low-yield emitters are presented below. 

Let us now turn to the role of interactions between 
dipoles in the ensemble, which is the main subject of this 
paper. Interactions play critical role in cooperative emis- 
sion since they introduce a disorder into system energy 
spectrum by causing random shifts of individual dipole 
frequencies;^^— In the conventional cooperative emis- 
sion, the main disorder effect is to split the narrow sub- 
radiant peak in the ensemble emission spectra. In the 
presence of metal nanostructure, radiation of subradiant 
states is expected to be quenched by much faster non- 
radiative losses in the metal. The crucial question is, 
however, whether interactions between closely spaced in- 
dividual dipoles can significantly alter the structure of 
collective eigenstates. In the remaining part of the pa- 
per, we present the results of our numerical simulations 
of cooperative emission fully incorporating both direct 
and plasmon-mediated interactions. 

We consider an ensemble of N molecular dyes attached 
to an Ag spherical particle with radius R — 20 nm via 
molecular linkers with approximately same length. The 
system is embedded in aqueous solution with dielectric 
constant eo = 1.77, and two types of dyes with quantum 
efficiencies q = 0.3 and q = 0.95 are used in the calcula- 
tions. A distinguishing feature of this system is a strong 
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FIG. 2: (Color online) Distribution of energy shifts for 20 
dipoles in C20 configuration around Ag NP at several dis- 
tances to its surface. 
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FIG. 3: (Color online) Distribution of decay rates for 20 
dipoles in C20 configuration around Ag NP at several dis- 
tances to its surface. 



effect of interactions on its geometryJ^ The flexible linker 
molecules hold the attached dyes with certain orientation 
of their dipole moments, so that repulsive inter-molecule 
interactions compel the dyes to form a spatially ordered 
structure on spherical surface. In our simulations, the 
dyes with normal dipole orientations were located at the 
sites of spherical lattice, specifically, fuUerenes C20, C32, 
C60, and C80, and, in some calculations, we included 
random deviations from the ideal lattice positions. The 
system eigenstates are found by numerical diagonaliza- 
tion of self-energy matrix, S^fc — Ajk — ^^jk, with its 
real and imaginary parts given by Eqs. (|M)) and (|^ . re- 
spectively. Calculations were carried at the SP energy of 
3.0 eV, the size-dependent Landau damping was incor- 
porated for all plasmon modes, and NP polarizabilities, 
Eq. ()13|) . with angular momenta up to / = 30, were calcu- 
lated using the experimental bulk Ag complex dielectric 
function. 

Figures [2] and [3] show distribution of real and imag- 
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FIG. 4: (Color online) Distribution of energy shifts for 20 
dipoles in C20 configuration around Ag NP at several average 
(with 10% fluctuations) distances to its surface. 



FIG. 5: (Color online) Distribution of decay rates for 20 
dipoles in C20 configuration around Ag NP at several average 
(with 10% fluctuations) distances to its surface. 



inary parts of complex eigenvalues of T,jk for N = 20 
molecules at the sites of C20 fullerene at three different 
molecule-surface distances. The system spectrum repre- 
sents several sets of degenerate eigenvalues indicating a 
high degree of lattice symmetry. For all distances, the 
histograms shovi^ a single eigenvalue vifith a large posi- 
tive energy shift (Fig. [2|), which corresponds to the di- 
rect dipole-dipole interaction between nearest-neighbor 
molecules. On the other hand, there are three degenerate 
eigenvalues with the largest decay rate, corresponding to 
predominantly superradiant states while the smaller de- 
cay rates are those of predominantly subradiant states 
(Fig. [3]). With increasing distance, the mixing between 
superradiant and subradiant states decreases and, for 
d — ZQ nm, decay rates of all but three eigenstates nearly 
vanish; note that in our approximation, pure subradiant 
states should have zero decay rate. Such behavior is due 
to diminishing contribution of higher-Z plasmons at larger 
distances [see Eqs. ( d^^ andlM])]. Importantly, direct in- 
teractions between close molecules result only in energy 
shift of subradiant states, without affecting superradiant 
states. We therefore conclude that dipole-dipole interac- 
tions do not destroy cooperative emission in plasmonic 
systems. 

This main conclusion remains unchanged when fluctu- 
ations (up to 10%) of molecules positions in radial direc- 
tion are included into simulations (Figs. S] and [S]). The 
spatial disorder lifts lattice symmetry, so that superradi- 
ant states now have different, however close, decay rates. 
Note that without interactions, i.e., when molecules an- 
gular positions are completely random, the spread of su- 
perradiant decay rates is considerably higher 

To elucidate the structure of collective states, we cal- 
culate the distance dependences of complex eigenvalues, 
Aj - fFj, for C20, C60, and C80 configurations of dyes, 
as shown in Figs. [Bl El and HI respectively. For C20 con- 
figuration (Fig. [5]) , there are five sets of eigenvalues with 
3, 4, 4, 7, and 1-fold degeneracies, in descending order 
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vs. distance for 20 dipoles in C20 configuration around Ag 
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sponds to three degenerate superradiant states, the dashed 
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of T J magnitudes. Down to the distance of c? = 5 nm, 
the largest decay rates, corresponding to three predom- 
inantly superradiant states, are well separated from the 
rest. The steep rise of Fj at small distances is due to 
increasing contribution of high-/ plasmon modes close to 
NP surface [see Eqs. and ^^]. The interplay be- 
tween various coupling mechanisms is especially revealing 
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FIG. 8: (Color online) Same as in Fig. [G] but for 80 dipoles 
in C80 configuration. 



when comparing the plots for Aj and T j (curves for same 
eigenvalue sets have similar patterns). By their d depen- 
dence, the eigenvalues fall into three main groups. The 
superradiant states have the largest decay rate Fj for 
all d and relatively small mainly positive frequency shift 
for d > i?/2; these states are dominated by plasmon- 
enhanced radiative coupling. The non-degenerate state 
with large positive energy shift and smallest decay rate is 
dominated by direct nearest-neighbor dipoles interaction; 
this state is least affected by the presence of NP and does 
not participate in the emission. The third group of states 
with mostly negative Aj and small Tj is dominated by 
nonradiative plasmon coupling. Closer to NP surface, 
the coupling becomes dominant due to high-^ plasmons 
and all states develop large decay rates and negative en- 
ergy shifts. Note that down to d > i?/4, the admix- 
ture between superradiant and subradiant modes is still 
relatively weak; below -R/4, the non-radiative coupling 
dominates the spectrum and the admixture is strong. 

For larger ensembles, the eigenstates have similar 
structure, as illustrated in Figs. [7] and [8] which show cal- 
culated eigenvalues for dipoles in C60 and C80 config- 
urations, respectively. Importantly, even with decreas- 
ing distance between the emitters in large ensembles, the 
dipole-dipole interactions still do not destroy coopera- 
tive emission. This can be understood from the follow- 
ing argument Mixing of superradiant and subradiant 
states takes place if the interactions between them are 
sufficiently strong. The latter requires that the electric 
field of a collective state is strongly inhomogeneous in 
space since, e.g., subradiant states couple only weakly 
to homogeneous field. On the other hand, such a field 



is comprised of individual fields of all the constituent 
dipoles so the resulting field's spatial fluctuations are 
weak if no two dipoles approach too close to each other, 
i.e., deviations of nearest- neighbor separations from their 
average, s = LN~^^'^, L being characteristic system 
size, are small. However, if deviations from s are large, 
i.e., two dipoles can be separated by a much closer dis- 
tance, s <^ LN^^/^, causing a strong spatial field fluctua- 
tion, then the eigenstates are no longer superradiant and 
subradiant states and cooperative emission is destroyed. 
This argument was confirmed numerically here by finding 
system eigenstates for both cases - dipoles on a spher- 
ical lattice with some fluctuations (see Figs. |l]and[5|), 
and a completely random angular distribution with fixed 
dipole-NP distance with no minimal separation between 
two dipoles (not shown). In the latter case, no superra- 
diant states were formed and the reason was traced to 
configurations with extremely close dipoles. Note, how- 
ever, that with both radial and angular distributions be- 
ing random, these are rare events. In the case of repulsive 
interactions between individual dipoles, considered here, 
deviations from the average dipole-dipole separation s 
are exceedingly small and cooperative emission survives 
the interactions. 

Another sharp contrast between plasmonic and pho- 
tonic Dicke effects is the fate of subradiant states. In the 
latter, the energy trapped in subradiant states is eventu- 
ally radiated, albeit with a much slower rate, resulting in 
sharp spectral features of emission spectrumj^i^ Instead, 
in plasmonic systems, the trapped energy is dissipated in 
the NP and only a small fraction of total energy leaves 
the system via superradoant states. Thus, the net effect 
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of plasmonic Dicke effect is to drastically reduce the emis- 
sion as compared to same number of individual dipoles. 
Remarkably, as the eigenvalues scale uniformly with N ^ 
the quantum efficiencies of superradiant states are nearly 
independent of the ensemble size, leading to the simple 
relation ([T]) that holds in the cooperative regime. 

This is illustrated in Figs. IHl and [TUl which show en- 
semble quantum efficiencies Qens [see Eq. p7p ] for two 
types of dyes with quantum yields q = 0.95 and q = 0.3, 
respectively. Two regimes can be clearly distinguished in 
the distance dependence of Qens'- it first shows a sharp 
rise with its slope proportional to N (non-cooperative 
regime) followed by a slower d dependence (cooperative 
regime). The crossover between two regimes takes place 
at d ~ 5 nm due to diminished high-Z plasmons con- 
tribution to nonradiative coupling for larger distances. 
In the cooperative regime, the precise behavior of Qens 
is affected by molecules' quantum yield. For high-g 
molecules, all Qens dependences collapse onto a single 
curve, Qens = 3(5, while for low-g molecules, Qens shows 



a weak dependence on TV. In both cases, this behavior 
can be easily understood from Eq. (|321)- Indeed, for de- 
generate superradiant eigenvalues we have = 1, and 
for large distances, as F — Fq, we obtain 

3 

^^"^ - l + 3(g-i-l)/iV' ^^^^ 

i.e., for large ensembles, the role of molecular quantum 
yield is diminished. 



V. CONCLUSIONS 

We studied here plasmon-mediated superradiance from 
an ensemble of dipoles near metal nanoparticle support- 
ing localized surface plasmon, thereby extending the 
Dicke effect to plasmonic systems. Our main conclu- 
sion is that the plasmonic Dicke effect is a robust phe- 
nomenon, more so than the usual photonic Dicke effect, 
because of a more efficient hybridization of individual 
dipoles via nanoparticle plasmon. We have established 
that hybridization takes place through two types of plas- 
monic coupling mechanisms - plasmon-enhanced radia- 
tive coupling and nonradiative plasmon coupling, the lat- 
ter having no analogue in the usual Dicke effect and caus- 
ing demise of cooperative emission at very close dipole- 
nanoparticle distances. 

While we considered a specific nanostructure - spheri- 
cal metal particle - the plasmon-mediated superradiance 
is a quite general phenomenon that should take place 
in any plasmonic system tuned into resonance with emit- 
ters, and Eq. ([T|) should apply provided that the usual cri- 
teria for cooperative emission are met. In fact, our theory 
remains unchanged for any nanostructure with spherical 
symmetry, for example metal nanoshells with dielectric 
core, upon simple replacement of NP polarizabilities in 
self-energy matrix Eq. (|12p and elsewhere with appropri- 
ate expressions. In fact, one expects that in nanoshells 
the nonradiatice losses would be smaller and so the plas- 
monic Dicke effect would be more robust than for solid 
nanoparticles. 

This work was supported in part by the NSF under 
Grant Nos. DMR-0906945 and HRD-0833178, and under 
the EPSCOR Program. 



Appendix 

Here we collect relevant some formulas for the elec- 
tric field Green dyadic in the presence of metal NP. The 
Green dyadic satisfies Maxwell equation 

V X V X G - fc2e(r)G = i, (45) 

where e(r) — e{uj)9{R — r) + eo9{r — R) is local dielectric 
function {0{x) is the step- function) . The Green dyadic 
can be split into free space and Mie-scattered parts. 
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G'^i,(r, r') = G°^(r, r') + G'^^(r, r'), where the free-space 
Green dyadic is 



G°.(r-r')=(^,.-^ 



)ff(r-r'), 



with 



47rr 

satisfying a scalar equation 

(A + fc2)g(r) = -5{r). 



(46) 



(47) 



(48) 



Consider first the free space part. Its near field expression 
can be obtained in the long wave approximation, i.e. by 
expanding in fcr ^ 1. In the first order, 



G° (r) 



1 li^aYu 



ik 



(49) 



In the far field limit, i.e., kr ^ 1 and kr' ^ 1, the free- 
space part can be expanded via Bessel functions. 



'ik\r — r 



-=ikY,Ji{kr')hi{kT)YUr)YL{r'). (50) 



Im 



47r|r — r'l 

which are approximated as 
[krj 



^ikr 



ji{kr') = 
yielding 



(2^ + 1).!' '^'(^^) = H)-V 



(51) 



G%{r,r') ^ (^S, 



ikr 

r 



1 

47r 



ikr' 



^yi,„(r)ri*„(r') 



After differentiation, the far field asymptotics takes the 
form 



A-TT 



(53) 



where we introduced xrm(r) = Vp[r'yim(f) andV';m(i") = 

v^[r-'-irz„(f)]. 

Now turn to the scattered part of the Green dyadic 
derived from solution of Mie problem for electromagnetic 
wave scattered on single sphere, 

G^,(r,r',k) = zkY,[<^iNL{^)Nrra{^') + 



Im 



contribution in Eg. ([54]) can be neglected as hi ^ li^^ 
The Mie coefficient a; has a form 

' tohi{pQ)[pii{p)]' - eji{p)[pohi{po)]' 

where pi — kiR, ki = ^^/el, and i — (e, eo)- For kR ^ 1, 
it becomes 



0,1 = — zs/a/Zc^'"''^, s; 



Z(2Z + 1)[(2;-1)!!]2' ^ ' 



where 



ai 



1 — isifc2i+iQ,j ■ 



(57) 



is NP multipolar polarizability that accounts for plasmon 
radiative decay, and 



_ p2;+i Ke-ep) 



ai = R 



(58) 



is the standard NP polarizability. The function Ni„i(r) 
is given by 



N,™(r) 



1 



k^l{l + l) 



V x[hf\kr)\.Yira{v)\, (59) 



where L = — i(r x V) is angular momentum operator. 
Using the following identity, 

V X [/i[i'(fcr)Lll,„(f)]=zrfc2 (fcr)l^„(r) 

+ iV[[krhf^'{kr) + h^^\kT)]Yi^{i)\, (60) 

prime standing for derivative, and expanding ft,|^-'(/cr) = 



(52) ii{kr) + ini{kr) in kr as 



h^Pikr) 



(kr)' (2Z-1)!! 
i- 



(2^ + 1)!! (fcr)'+i 



(61) 



we obtain 

^hn{kr) 

where 



ky/si{2l+l) 



V[^iikr)YU^)], (62) 



(pi{kr) 



[kr) 



isi{kr)K 



(63) 



Thus, for kr and kr' ^ 1, the scattered part of the 
Green dyadic has the form 



biMl'^{v)MUv')\, (54) G^.(r,r',k) « z/c^[a,iV,';„(/cr)iVL(A:r')] 



(64) 



Im 



where the first and second terms are electric and mag- 
netic contributions and ai and 6; are the Mie coefficients. 
In the long wave approximation, kR <^ 1, the magnetic 



Y.^^,[vi{kr)YU^)]K[Mkr')Y;^{r')]- 



Im 



12 



This expression can be further simphfied by substituting 

ai = ai +isik^^+^\ai\'^, where ai = a,|l - isiP^+'^ai\~^ , 
and keeping the first two powers of k 

Im 

E [5iK™(r)xr™(r')+x5'™(r)^r™(r')] 



m— — 1 



ai|>f,„(r)Vr™(r') 



(65) 



which, after adding the free-space part of the Green 
dyadic and neglecting plasmon radiative decay, leads to 

Eq. dm). 



For fcr > 1, with help of Eqs. ([5Tj), ([591), and (gOl), we 
easily obtain 

ikr 

Nz„(r) = - , , Vyzn.(r), (66) 

A:^/(« + 1) 

and combining this expression with Eq. (j62p . we obtain 
the far field Green dyadic (i.e., fcr ^ 1 and kr' <C 1) 



G^.(r,r') = -^ 5] e^'='-[V^ri„(f)]^r™(r'), (67) 



where we set I = \. 
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